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rn ■ ABSTRACT 

In Spitzer observations of Tauri stars and their dislcs, PAH features are detected in less than 10% 
I of the objects, although the stellar photosphere is sufficiently hot to excite PAHs. To explain the 

' deficiency, we discuss PAH destruction by photons assuming that the star has beside its photospheric 

(-H ' emission also a FUV, an EUV and an X-ray component with fractional luminosity of 1%, 0.1% and 

C^' 0.025%, respectively. As PAH destruction process we consider unimolecular dissociation and present a 

Q ■ simplified scheme to estimate the location from the star where the molecules become photo-stable. We 

, find that soft photons with energies below ~ 20 eV dissociate PAHs only up to short distances from 

, the star (r < lAU); whereas dissociation by hard photons (EUV and X-ray) is so efficient that it 

would destroy all PAHs (from regions in the disk where they could be excited). As a possible path for 
PAH survival we suggest turbulent motions in the disk. They can replenish PAHs or remove them from 
' the reach of hard photons. For standard disk models, where the surface density changes like r and 

' the mid plane temperature like r""'^, the critical vertical velocity for PAH survival is proportional to 

' y.-3/4 and equals ~5 m/s at 10 AU which is in the range of expected velocities in the surface layer. The 

uncertainty in the parameters is large enough to explain both detection and non-detection of PAHs. 
'_ ' , Our approximate treatment also takes into account the presence of gas which, at the top of the disk, 

G^ . 

, is ionized and at lower levels neutral. 

Key words, dust, extinction - planetary systems: protoplanetary disks - infrared: stars - X-rays: stars 
^ ■ - X-rays: ISM 

1. Introduction 

Infrared emission bands of PAHs can be used as a probe of the UV environment. They are commonly 
seen in the ISM, but also in young stellar objects such as Herbig Ae/Be stars (Waelkens et al. 1996, 
Siebenmorgen et al. 2000, Meeus et al. 2001, Peeters et al. 2002, van Boekel et al. 2004). The ob- 
served emission can be explained in models of an irradiated disk (Habart et al. 2004, Visser et al., 2007, 
Dullemond et al. 2007a). 

ISO also looked at a few of the much fainter T Tauri stars but without a clear PAH detection 
(Siebenmorgen et al. 2000). In the Evans ct al. (2003) legacy program which employs the more sensitive 
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Spitzer Space Telescope (SST), 3 out of 38 T Tauri stars show PAH features (Geers et al. 2006). This 
corresponds to a detection rate of only 8% in contrast to almost 60% in Herbig Ae/Be stars (Acke & van 
den Ancker 2004). Similarly low rates for T Tauri stars are found by Furlan et al. (2006) who present 111 
SST spectra in the Taurus- Auriga star forming region and speculate that the absence of PAH resonances 
is due to the much weaker UV field compared to Herbig Ae/Be stars. Geers et al. (2009), on the other 
hand, argue that the PAHs are simply under-abundant relative to the ISM. They also find that variations 
of the disk geometry, such as flaring or gaps, have only a small effect on the strength of the PAH bands. 
Clearing out gas and dust by planet formation inside the disk could effectively remove PAHs. Indeed, 
inner gaps in disks are observed at radii between 40-60 AU and at wavelengths between 20-1000/^m where 
the emission is dominated by large grains (Brown et al. 2008, Geers et al. 2007b). However, in cases where 
PAH emission is resolved, it is extended up to 15-60 AU, without sub-structure and inside the inner gap 
region (Geers et al. 2007a). The spatial extent of the PAH emission is also similar for T Tauri and Herbig 
Ae/Be stars. We therefore suggest that PAH removal by radiative destruction is dominant. 

Present radiative transfer models of the PAH emission from dusty disks consider only the stellar 
radiation field and no additional EUV or X-ray component (Habart et al. 2004, Geers et al. 2006, Visser 
et al. 2007, DuUemond et al. 2007). Their hard photons could, according to laboratory experiments (Ruhl 
et al. 1989, Leach et al. 1989a, b, Jochims et al. 1994) and theory (Omont 1986, Tielens 2005, Rapacioli 
et al. 2006, Micelotta et al. 2009), destroy PAHs. We discuss below their impact on the PAH abundance 
in the disks of T Tauri stars. 

2. Radiation components of T Tauri stars 

Our T Tauri model star has a total luminosity = 2Lq. Its radiation consists of a photospheric, a 
FUV, an EUV and an X-ray component. Their parameters are listed in Table [T] and are very similar to 
those proposed by Gorti & HoUenbach (2008). The total spectrum is displayed in Fig[T] We point out 
that the FUV and EUV radiation are observationally poorly constrained. 

The photosphere supplies most of the luminosity whereas the FUV, EUV and X-ray radiation are 
much weaker and believed to originate from accretion onto the star and from chromospheric and coronal 
activity. The photosphere, the FUV and EUV component are approximated by blackbodies. We assume 
4000 K for the photosphere and, following Stabler et al. (1980) and Calvet & GuUbring (1998), 15000 K 
for the FUV (pre-shock) and '^3xlO^K for the EUV emission (post-shock region). 

The strength of the FUV and EUV radiation is determined by the accretion luminosity which we 
approximate by Lace — GM^M / R^. If = 2R0 and Af* — IMq are the radius and mass of the 
star, an accretion rate M = 10^^ Mq yr^^ (Akeson et al. 2005) yields iacc ^ 0.01 L*. Higher values 
[Lacc/L^, > 0.1), but with a large spread, are derived by MuzeroUe et al. (1998, 2003) from hydrogen 
emission lines. However, as we show in section[51 such stronger fluxes have little influence on the stability 
analysis of PAHs. 

Preibisch et al. (2006) establish from Chandra observations (0.5-8keV) a relation between the X- 
ray luminosity, Lx^ and the total luminosity confirming the ROSAT results of Sterzik & Schmitt 
(1997). The ratio L^/L^ is similar in rapidly rotating main-sequence stars and non-accreting T Tauri 
stars (~10^'^), but systematically lower by a factor ~4 in accreting T Tauri stars (Preibisch et al. 2006). 
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Fig. 1. The spectral energy distribution of our T Tauri model star at 1 AU without foreground extinction 
(Eq. [1]) . The absolute luminosities of the components are given in Table [TJ 

Table 1. The four radiation components of our T Tauri model star. 
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(1) fractional luminosity, 

(2) spectral shape (BB — blackbody), 

(3) mean energy of emitted photons, hV^ 

(4) approximate flux at 10 AU, Fio. 



Interestingly, in Herbig Ae/Be stars L^/L^, is much smaller (^10^^ . . . 10^^, Stelzer et al. 2006) with 
values comparable to the Sun. The solar X-ray luminosity in the 0.1-2.4keV ROSAT passband lies during 
a solar cycle in the range 10^^'* < L^/L^, < 10^^'^ and is typical for G stars (Judge et al. 2003). 

X-ray fluxes are generally variable on timescales of hours to weeks and weaken during the evolution of 
the T Tauri star. For example, half of the sources in the Taurus molecular cloud detected by XMM/Newton 
(0.3 - 7.8 keV) show variations, more at hard (> 0.5 keV) than at soft energies, and a quarter of them 
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display flares (Stelzer et al. 2007), about once a week and lasting for a few hours. In a strong flare, more 
than lO'^^ erg are emitted and can reach 1% of the total luminosity. We assume up to 2keV a power 
law spectrum oc (Giidel et al. 2007) and neglect harder radiation because the emission then steeply 
declines (cx v^^). 

Let Li be the frequency-integrated luminosity of the radiation component i (see Table [Ij and L^^^ its 
spectral luminosity such that Li — J Li^^, dv. Dropping for convenience the index i, the flux (of component 
i ) at a distance r is 

= W 

where we included a screening factor e~'^" to account for foreground absorption (by dust and gas). If 
denotes the absorption cross section per carbon atom, a PAH of Nc carbon atoms absorbs in one second 
(from component i) 

N^^N.j?^,. (2) 

photons of total energy 

Eahs = y F^i^v dv (3) 
The inverse of is the average time between two absorption events, 

iabs = N-^ (4) 
The mean photon energy equals 

J hu 

3. Cross sections 

As the light from the star enters the disk, it is attenuated by gas and dust. The absorption cross sec- 
tion of gas depends on the ionization stage of the atoms which is determined by the balance between 
recombination and photo-ionization. By far the most important atoms are, of course, hydrogen and he- 
lium with ionization potentials of 13.6 eV and 24.6 eV, respectively. Because the recombination rate is 
proportional to the square of the gas density which is high in the disk (section 5), the gas is ionized 
only in a thin surface layer {Ay < 0.001 mag, section 2.6 of Gorti & HoUenbach 2008). We use atomic 
cross sections of Morrison & McCommon (1982) and Balucinska-Church & McCommon (1992) and solar 
element abundances. 

The dust cross sections are taken from the model of Kriigel ( 2006) which describes standard dust. For 
X-rays, the absorption efficiency calculated from Mie theory must be corrected downwards. Hard photons 
can eject electrons from the grain and as these carry away kinetic energy, only part of the photon energy 
is deposited in the dust particle. The threshold, Et, above which such a correction is necessary depends 
on the grain size; details are given in Dwek & Smith (1996). For a 10 A graphite particle, Et ~ 100 eV 
and the reduction factor is roughly proportional to v^^. 
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Fig. 2. The mass extinction coefficient per gram disk material when the gas is neutral (Morrison & 
McCommon 1983); the gas-to-dust mass ratio equals 130. 

The absorption coefficient of dust, K^^x, and of neutral gas plus dust, K\ = Kg^s.x + -K'd.A, both 
per gram of disk material, are plotted in Fig. [2] for a dust-to-gas mass ratio of 1:130. Note that at the 
ionization threshold of hydrogen, ifgas is almost 10^ times greater than K^- 

With respect to the absorption cross section of PAHs, we assume = 7 x 10^^^ cm^ per carbon atom 
when hv < 13.6 eV and scale at higher energies to follow the values of a graphite sphere of 10 A radius 
(Dwek & Smith 1996). The maximum wavelength (in A) for PAH excitation is Amax.PAH — 1630-|-370\/]Vc 
(Schutte et al. 1993) resulting in a minimum photon energy of 2.3 eV for a PAH with Nc = 100 carbon 
atoms. 

3.1. PAH emission 

As the PAHs are transiently heated, their excitation is usually treated statistically. Following 
Guhathakurta & Drainc (1989), let P{T)dT be the probability of finding in a large ensemble of PAHs 
in a steady state an arbitrary PAH in the temperature interval [T, T + dT] . The temperature distribu- 
tion function P{T) is calculated in this method from a transition matrix {Afi). If denotes the PAH 
absorption cross section, the matrix element Afi referring to dust heating from an initial enthalpy bin 
centered at Ui to a final one centered at Uf and of width is, for a mono-chromatic flux, equal to 

^,.4^^ "'"'-''•-'"'l^^^''' (6) 
[ : else 

Examples of P{T) are displayed in Fig. [3] for mono-chromatic fluxes which cover almost the en- 
tire range encountered anywhere in the disk in terms of intensity, the flux ranges from = 10 to 
10*" ergs~^cm~^, as well as hardness, the photon energy is between hv — 3.8 eV and IkeV. The approxi- 
mate unattenuated fluxes of the four radiation components at a distance of 10 AU are listed in Table [TJ 
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Note that when F — n-^hv is constant, the number of photons n-^ decreases as the photon energy hv goes 
up. The power absorbed by one PAH is ahiiost independent of the photon energy as long as hv < 100 eV 
(see Fig. By and large, when hv is fixed and F increases, the curves in Fig. [3] narrow and move to 
the right towards higher temperatures. When the radiation field is weak {F — 10ergs~^cm~'^) and the 
photons are soft {hv = 3.8 eV), the PAH absorbs about one photon per day and there is plenty of time 
to cool down. In this case, the PAH virtually never exceeds the sublimation temperature Tg, (at which 
solid carbon gasify). For a wide pressure range (10~^ — 10~^dyn cm~^), Tg ^ 2000K for graphite (CRC 
Handbook of Chemistry & Physics 2005, Salpeter et al. 1977). In Fig. [3] we highlight the area where the 
temperature is above Tg. For a PAH exposed to a strong radiation field {F = 10'' ergs~^cm~^), about 
a dozen of soft photons are absorbed within a cooling time and the temperature distribution function 
becomes very narrow around 2000 K and sublimation is likely. For hard photons {hv > 50 eV) , the PAH 
undergoes, independent of the strength of the radiation field, extreme temperature excursions. This indi- 
cates that PAH become photo-unstable either by absorption of a single hard photon, or by soft photons 
if there are many of them. 

4. PAH destruction 

The abundance of PAHs is determined by the competition between formation and destruction processes 
under the specific environmental conditions. Underlying processes are discussed, for example, by Omont 
(1986), Voit (1992), or recently by Micelotta et al. (2009a). Here we only consider PAH destruction 
by photons and generally assume that PAH formation is negligible. After photon absorption, a highly 
vibrationally excited PAH may relax through emission of IR photons or, if sufficiently excited, lose 
atoms. The latter process is called unimolecular dissociation and is discussed for interstellar PAHs by 
Allamandola (1989), Leger et al. (1989), Le Page et al. (2003), RapacioU et al. (2006), and Micelotta et 
al. (2009b). Laboratory studies of PAH dissociation which can be applied to astrophysical situations are 
rare (Jochims et al. 1994). The photo-chemistry of PAHs is reviewed by Tielens (2005, 2008). 

4.1. Procedure 

In the disks of T Tau stars, the PAH abundance depends obviously on place and on time as the disk 
evolves. There is no general solution to the problem and to extract numbers, we have to radically simplify 
it. We wish to find some estimate of the location where PAHs become stable against photo-destruction. 
To derive a procedure, we recall that although after absorption of an energetic photon its energy is 
immediately distributed over all available vibrational modes (Allamandola et al., 1989), the excitation 
of a particular atom fiuctuates and occasionally it is pushed into the continuum and leaves the PAH. 
Quantitatively, the unimolecular dissociation can be written in Arrhenius form. In a classical description, 
an atom of critical (Arrhenius) energy Eq detaches from a PAH of peak temperature Tp if the dissociation 
time 

idis - v-'e^^/'^^" (7) 

is shorter than the cooling time icooi- A characteristic value for the vibrational frequency is vq ^ 
10^"^ s~^. The "atom", which may besides H or C also be an atomic group like C2H2, needs the time idis 
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temperature [K] 

Fig. 3. The temperature distribution P{T) of a PAH with 100 C atoms exposed to mono-chromatic 
radiation with hi' = 3.8, 13.6, 50, 100, 300eV and 1 keV. This set includes the mean photon energies of 
the four radiation components of the T Tauri star. The fluxes range from (top to bottom) F = 10 to 10^ 
erg cm~^. Shaded area marks temperatures above the sublimation temperature of graphite. 



to overcome the critical internal barrier, Eq, which is similar but not identical to the chemical binding 
energy. Micelotta et al. (2009b) quote Eq of 3.2eV for H-loss, 4.2eV for C2H2, 7.5eV for pure C loss and 
9.5eV for C2. For the ISM they find Eq = 4.6 eV and a somewhat larger value for PDR. The inverse of 
the dissociation time is the probability that a certain atom leaves the PAH per unit time. 
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Fig. 4. The power W ^ Nc k F absorbed by one PAH, with number of C atoms Nc as indicated, in a 
mono-chromatic flux F — 10* erg s^^ cm^^ as a function of photon energy = hv. 

The exponential term e^°l^'^v in Eq.© increases very rapidly as T falls and meaningful values (i.e. not 
too large ones) of idis are obtained only if Tp > 1500 K. Atoms will only detach when tdis < ^cooi- As the 
cooling time at these temperatures is for astrophysical applications of order 1 s, independent of the PAH 
size, the dissociation criterion reads 



Id 



It leads to a minimum temperature for destruction 

^0 



(8) 



(9) 



fc In vq 

Assuming £^0 ^ 5eV, one gets Tdis — 2000 K. At this high temperature, the internal energy of a PAH 
is reasonably well approximated by SiVcfcTdis when also taking the presence of H-modes into account. 
The minimum temperature Tdis is related to a minimum energy input IS.E. The PAH is therefore unstable 
to photons with 



AS > hv^ = 37VcA:Tdis = 



In i/Q ' 



or when the number of carbon atoms 
2b.E 



Nr < 



[eV] 



(10) 



(11) 



Micelotta et al. (2009b) find that a PAH with Nc — 50 requires an internal energy of about AE = 24 eV 
to dissociate, which agrees with the above estimate (Eg. [TT|) . 
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Fig. 5. Of each radiation component, ~90% is absorbed in what we caU the extinction layer. The optical 
depth from its bottom to the star is one and in vertical direction equal to the grazing angle ckg. The 
height of its lower boundary, zq, declines with radius, but its geometrical thickness is rather constant 
{£ ~ 0.5 H, see FiglHl). Vertical motions may replenish PAHs from below. 

The minimum energy input required for dissociation can either be delivered by absorption of many 
soft photons, with a total energy E'abs > Ai? (Eq.[3]), or ii) by a single hard photon, with energy hi^ > AE. 
If a photon heats the PAH to a peak temperature much above T^is, more than one atom will detach. The 
first expulsion occurs momentarily (tdis <C Is). It consumes the energy Eg plus some kinetic energy E^kin 
for the liberated atom. The new PAH temperature follows from 



hi^^ Eo- Ei,in = 3{Nc~l)kT 

This happens x times until T has dropped to Tdis, 

hv -x{Eq + Skin) = 3(7Vc - x)kTdis 



(12) 
(13) 



With i?kin ^ 0.5 eV, we estimate that the total number of freed atoms is 

_ hi' ~ SNckTdis hv Nc , , 

Eo + Ei^in - SkTdis ~ 5 [eV] " To ^ ^ 

For Nc — 100, an average EUV photon ejects nine atoms and an X-ray photon destroys the whole PAH 
(column (5) in Table [2]). 

4.2. Disruption by Coulomb forces 

For completeness, we also mention disruption of PAHs by Coulomb forces. Double or multiple ionization 
of a PAH loosens the binding of the peripheral H atoms as well as of the skeleton of carbon atoms. The 
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ejection of K-shell electrons by X-ray photons {hv > 284 eV) in combination with Auger electrons will 
amplify the process. Coulomb explosion is relevant mainly for small PAHs and neglected here. 



5. Conditions for PAH survival 

According to Eq. pI7|) . PAHs are destroyed if the source emits photons of energy hv > O.INcEq, irre- 
spective of the distance to the star or its luminosity. For Nc ~ 100, the critical photon energy is only 
50 eV (Eq fTU)) . As T Tauri stars (or their jets) also radiate at X-rays and in the EUV, the surface of the 
disk should be devoid of PAHs unless a) the period over which hard photons are emitted is too short 
to destroy all PAHs; b) the PAHs are by vertical motions removed from the hard radiation before they 
are destroyed and there is an influx of PAHs from below; c) PAH destruction is compensated by PAH 
formation in the surface layer. The last effect should, in a hard photon environment where PAHs and 
carbon atoms are ionized, be prohibited by Coulomb repulsion (Voit, 1992). 

5.1. Destruction time 

The above PAH survival condition under a) can easily be dismissed. To estimate the time for PAH 
removal, ircm, by the radiation component i, we note that most of the radiation is absorbed on the disk 
surface in a sheet of vertical optical depth t± equal to the grazing angle of the incident light. We call 
this sheet the extinction layer (of radiation component i) and denote its geometrical thickness £i (see 
Fig[5]). To first order, the PAHs in the extinction layer receive the stellar flux of Eq.((T]) with = 1. If 
the instability criterion of Eg. pTI)) is fulfilled, iicm follows from xN^tj-cm = Nc, where x is from Eq. p^ . 
therefore 



With iabs from Tabled one sees that even at 100 AU, ircm is short compared to the duration of the 
T Tauri phase (-10^ yr, Bertout et al. 2007, Cieza et al. 2007). 

5.2. Exposure time and vertical mixing 

Next we consider the possibility that vertical motions in the disk lead to a continuous exchange between 
matter in the extinction layers, where almost all photons are absorbed and PAHs destroyed, and the 
layers below where PAHs are shielded and damaged ones possibly rebuilt (Fig[5]). We assume that gas 
and dust are perfectly mixed in a mass ratio 130:1. 

In a Keplerian disk that is isothermal in z-direction and in hydrostatic equilibrium, the gas density 
changes like 



'rem 




(15) 




(16) 



Here S(r) is the surface density at radius r which is assumed to follow a power law. 



oo 




7 



(17) 







R. Siebenmorgen and E. Kriigel: PAH destruction and survival in T Tauri disks 



11 




Fig. 6. The height, zq, of the bottom of the extinction layer and its thickness £ for the four radiation 
components (see Figl5]and Eg. ([211) . ([22]) ) for a grazing angle Ug = 3°. Due to the high gas densities, the 
EUV extinction layer is practically coincident with the ionization front and lies above the extinction layer 
of the other components. This implies that EUV radiation is absorbed first whereas the other components 
penetrate deeper. Neglecting ionization by the FUV component, the gas below the ZQ-line of the EUV 
radiation (top) is neutral. 

and 

H{r) = y/kTr^/GA'Um (18) 

is the scale height, Af* ~ 1 Mq the stellar mass and m the mass of a gas molecule. For the surface 
density, reasonable numbers are 7=1 and S = 200 gcm^^ (Hartmann et al. 1998, Kitamura et al. 2002, 
DuUemond et al. 2002, Rafikov & CoUe 2006, Gorti & HoUenbach 2008), although the various estimates 
show considerable scatter. 

For the radial variation of the gas temperature in the opaque mid plane, T(r), we also adopt a power 
law, 

-0 

(19) 

The mid plane is roughly isothermal in z because the optical depth is high and the net flux zero. It is 
much colder than the extinction layers because it is not exposed to direct stellar heating. The radiative 
transfer in the disk, including the energy equation, can be solved to any desired accuracy even when the 
disk is very opaque (see section 11.3.2 of Kriigel 2006). As long as the dust in the mid plane is optically 
thick to its own emission, the results for T{r) can be well approximated by putting in Eg.p^ (3 = 0.5 
and To ^ 130 K (as also suggested by DuUemond et al. 2007b or Chiang & Goldreich 1997). 

Each extinction layer extends vertically from some value zq upwards to infinity (Figl5]). We give it a 
finite thickness £ by demanding that, say, 90% of the photons are absorbed between zq and zo + L If v± 
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Fig. 7. The critical vertical velocity for PAH survival «„• after Eq. (|23|) with respect to X-rays as a 
function of distance (full line). We identify the vertical velocity v± with the turbulent velocity Wf PAHs 
survive when v± — vt > Vcr, else they are destroyed (shaded areas) by expulsion of atoms. The hatched 
strip between the dashed lines shows the range where v± is between 0.001 Cg and 0.01 Cg (where Cg is the 
sound velocity). 



denotes the typical vertical velocity, for example, as a result of turbulence, PAHs are exposed to radiation 
for a time 

ioxp = — (20) 

This is also the mean residence time of a PAH in the extinction layer. For PAHs to survive, <oxp must 
be smaller than ij-cm Eq.ljTS]). The height zq follows from 

/>oo 

K p{z)dz^a^ (21) 

J Zo 

and £ may be estimated from the condition that only 10% of the flux is absorbed above zq + i, 



K I p{z)dz = 0.1 Qfg (22) 

K is the mass absorption coefficient of gas and dust at the characteristic frequency of the particular 
radiation component (see Table[2]). Because p{z) changes rapidly, is rather insensitive both to Ofg as well 
as K . For a^/K = 10~^ . . . 10~^, one obtains zq = 2.6H . . . 5.7H. So in the V band, where absorption 
is only by dust {K ~ 200 cm^ g~^) and for a grazing angle ag — 3°, one gets zq ~ {4...5)H and 
fi — 1/ H 0.5. The height zo where an extinction layer begins and its thickness are shown in FiglHl ^ 
is for all radiation components very similar (£; ~ H/2, i = 1, . . . , 4). EUV photons are absorbed highest 
up, their extinction layer lies about one scale height above the others. X-rays penetrate slightly deeper 
than photospheric or FUV photons. 
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Table 2. Quantities relevant to PAH survival. 
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hV 
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K/Ky 
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m/s 


photosphere 


2.7 


700 


< 1 


1.2 




4.4 


0.5 


5 s 




FUV 


4.2 


700 


< 1 


1.5 




4.5 


0.5 


120 s 




EUV 


97 


36 


~0.5 


110 


9 


5.3 


0.4 


9 days 


2900 


X-ray 


1100 


0.17 


< 1 


0.55 


100 


4.2 


0.5 


290 yr 


3.4 



(1) The mean energy of destructive photons, hV, from Eq.((5]) at the bottom of the extinction layer (see FigEJ; 
for photospheric and FUV photons we put i/c ~ and for EUV and X-rays we integrate for v > i/c (Eg. I10|l: 

(2) absorption cross section, k, per C atom at frequency 77; 

(3) approximate mean degree of ionization of the gas in the extinction layer; 

(4) extinction cross section of gas and dust at frequency V normalized to Kv = 200 cm^ g~^; 

(5) number of expelled atoms, x, per absorption event from Ea. p4|l : 

(6) altitude of the bottom of the extinction layer in units of the scale height H; 

(7) thickness of extinction layer; 

(8) mean time tabs from Eq.Q in which one photon is absorbed; 

(9) critical vertical velocity for PAH survival. 
Note: zo,i,tahs and Va refer to r = 10 AU. 



From Eq.Q, and (|20|) . one finds that vertical motions safeguard PAHs against destruction if 

with Vc from Eq.([lOl). When the scale height H is given by Eq.lUHl), oc r^^^+^^Z^ = rpj^g 
values of Wcr, at a distance of 10 AU, are listed in Table [5] together with other quantities relevant to PAH 
survival. 

One expects the disk also to be turbulent. Turbulence may be driven by various processes such as 
shear flows in the disk (Lin & Bodenheimer 1982), magneto-rotational instabilities (Balbus & Harley 
1991) or velocity discontinuities at places where infalling matter (Cassen & Mossman 1981) or outflows 
(Elmegreen 1978) strike the disk surface. Until now it is not clear which type of turbulence dominates. 
We assume that for the size of the largest Eddies, £cd, the average turbulent velocity, Vt, grows linearly 
with the sound speed Cg. Various hydrodynamical 3-dimensional calculations (Cabot 1996, Boss 2004, 
Johansen & Klahr 2005, Fromang & Papaloizou 2006) support this view. The favoured parametrisation 
is Vt = a'^Cg with q = 0.5. This choice has consequences on the Eddy scale t^d = a^~'^H and the turn over 



time ted = icd/vt — ^"^/f^K, with Keppler frequency Qk = y/GALjW (e.g. DuUemond & Domink 
2004). Weidenschilling & Cuzzi (1993) use g = 1 so that the Eddy scale is about the pressure scale height, 
€cd = H , and larger than the thickness of the extinction layer, ^ed > ^ ^ H/2. Estimated values for a 
are in the range from 0.0001 up to 0.1 (DuUemond & Dominik 2004, Schrapler & Hcnning 2004, Youdin 
& Lithwick 2007). Taking a — 0.01 the Eddy scale is 10 times smaller for q — 0.5 than for q — I, and 
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Fig. 8. The energy E^hs (Eql3]), which is absorbed by a PAH of Nc ~ 100 carbon atoms, as a function of 
distance from the star. The PAH is exposed to the photospheric and FUV radiation component described 
in Table [H The dashed hne refers to the top of the extinction layer (r — 0) and the full line to its bottom 
(r = 1). For Arrhenius energy of i?o = 5, 7 and 10 eV, the minimum energy input AE (EqllOp for PAH 
dissociation is indicated by the dotted lines. 



in addition larger turbulent velocities are obtained with vt = \/aCs, supporting a faster transport of the 
PAH. Identifying v± in Eq. (|23p with the turbulent velocity vt and assuming a temperature dependence 
as in Eq. fTO]) . we plot in FiglT] the vertical velocity v± as a function of radius. The figure also shows 
the critical velocity for PAH survival, Vcr, with respect to X-rays. Note that the critical velocity is for 
the X-ray radiation component insensitive to the particular choice of Eq and i^c (Eq. [TOl 1^5]) . When 
v± = 0.01 Cg, PAHs can survive at distances r > 10 AU; when v± is considerably smaller than 0.01 Cg, 
they cannot. 

Critical velocities for PAH survival are much higher for EUV than for X-ray photons (Tabled]). EUV 
radiation will therefore always destroy PAHs but, as depicted in Fig El the EUV extinction layer is the 
topmost and below it, PAHs may survive and be excited. 

If PAHs are removed from the extinction layer before they are destroyed, they must, in order to be 
detected, at the same rate be injected from below. Therefore, the critical velocity can alternatively be 
expressed through 




'cr 



(24) 



which leads to similar values. We note that the mass reservoir below the extinction layer is sufficient to 
sustain, over the lifetime of the disk tiifc, the required mass influx pvcr- 
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5.3. PAH dissociation by soft versus tiard ptiotons 

The energy i?abs (Eql21), which is absorbed by a PAH of Nc — 100 carbon atoms, is shown in Fig. [5] 
as a function of distance from the star. The PAH is exposed to the photosphcric and FUV radiation 
component described in Table [1] and results are shown for the top (r = 0) and the bottom (r = 1) 
of the extinction layer. The minimum energy input Ai? (Ea llOp for PAH dissociation depends on the 
choice of the Arrhcnius energy Eq and is indicated for Eq ~ and 10 eV, respectively. In this picture, 
for i?o = 5 eV and at the bottom of the extinction layer, PAHs are dissociated by soft photons up to 
1 AU. For X~rays, however, we find that PAH destruction occurs typically at distances up to ^ 10 AU or 
even larger (Fig[71). Dissociation of PAH acts for soft (photosphcric and FUV) photons on much shorter 
distances than for hard photons (X-ray component). 

6. Large grains 

Observations of T Tauri stars at millimeter wavelengths (Testi et al. 2003, Lommen et al. 2007) and in 
the mid infrared (van Boekel et al. 2003, Przygodda et al., 2003, Kessler-Silacci et al. 2006, Bouwman et 
al. 2008, Watson et al. 2009) suggest that grains in T Tauri disks are at least 10 times larger than those 
in the ISM. As such large grains may also be present in the top disk layer we estimate how this would 
affect the stability analysis of PAHs. We first note that in case of homogeneous mixing an increase in 
particle size would not alter the dust-to-gas mass ratio. 

Should the grains be much larger than the wavelength, the absorption coefficient per gram of dust, 
Kd,\i would decrease roughly like one over grain radius whereas the ratio K^.x/ K^y would still roughly 
be given by the values in Table [51 For hard X-rays, on the other hand, K^ x is not sensitive to grain size. 

Therefore, if disk grains are on average ten times bigger and thus 10'' more massive than interstellar 
ones, we expect that the height zq to which the stellar radiation components can penetrate (see Fig[5]) 
stays the same for X-rays, but also for EUV radiation because EUV absorption is due to gas, not dust. 
However, optical and FUV photons will reach farther down, about half a scale height, so that there 
may be a thin disk layer {^H/A) where PAHs are shielded from X-rays and EUV photo-destruction and 
excited by optical or FUV radiation. 

7. Conclusion 

In the search for an explanation why most T Tauri stars do not exhibit PAH features and only a few 
do, we investigate which processes can remove PAHs from the surface layer of T Tauri disks and under 
which conditions they should be present. Clearing of PAH through interaction with planets seems not an 
efficient process (Geers et al 2007b) and we show that PAH under-abundance can be caused by radiative 
destruction. We use a fiducial model for the photon emission of the T Tauri star that includes, beside 
the photosphere, FUV and EUV radiation and an X-ray component. 

1. We introduce for each stellar radiation component the notion of extinction layer as the place where 
^^90% of the photons are absorbed. EUV photons are mainly absorbed by gas, X-rays by gas and dust 
alike and the photosphcric and FUV component arc only attenuated by dust. The extinction layer of 
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all four components have a similar geometrical thickness and their bottom is at similar elevation zq, 
except for the EUV extinction layer which lies higher up (FigE]) . 

2. PAH may be radiatively destroyed, by unimolecular dissociation, where one or several atoms are 
expelled after photon absorption. 

3. Destruction by the photospheric and FUV radiation component (soft photons), increases with the 
strength of the radiation field and is very efficient below f AU. 

4. Hard photons can dissociate PAHs at all distances and their efficiency grows with the hardness of 
the photons. Without some counter process, all PAHs (in layers where they can be excited) would be 
destroyed within a time short compared to the lifetime of the disk. 

5. Although grains in the disk surface are presumably larger than interstellar ones, the stability analysis 
of PAHs would not change significantly. 

6. Therefore, in disks where PAHs are detected, there must be some survival channel. Because creation 
of PAHs in the extinction layer is too slow to compete with PAH destruction (Voit f 992), we suggest 
vertical mixing as a result of turbulence. It can replenish PAHs or remove them from the reach of 
hard photons. 

7. For standard disk models, the minimum velocity for PAH survival is proportional to r~^/^ and equals 
~5m/s at 10 AU. If turbulent velocities are proportional to the sound speed a velocity Ut > 5 m/s 
would imply Vt/cs <; 0.01 as PAH survival condition. Theoretical predictions for this ratio have a large 
spread but in accordance with the observational fact that PAH features are usually absent it seems 
that generally the condition is not fulfilled. 

8. A higher PAH detection rate is found in Herbig Ae/Be stars. In our picture this is explained as their 
destructive hard radiation component is relatively weak {Lx/L^, < 10^^, Preibisch et al. 2006) and 
because the intensity of the PAH emission from large distance from the star is larger given their 
higher optical luminosities. 

Acknowledgements. We thank the second anonymous referee for constructive comments. 
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ABSTRACT 

In Spitzer observations of Tauri stars and their disks, PAH features are detected in less than 10% 
of the objects, although the stellar photosphere is sufficiently hot to excite PAHs. To explain the 
deficiency, we discuss PAH destruction by photons assuming that the star has beside its photospheric 
emission also a FUV, an EUV and an X-ray component with fractional luminosity of 1%, 0.1% and 
0.025%, respectively. As PAH destruction process we consider unimolecular dissociation and present a 
simplified scheme to estimate the location from the star where the molecules become photo-stable. We 
find that soft photons with energies below ~ 20 eV dissociate PAHs only up to short distances from 
the star (r < lAU); whereas dissociation by hard photons (EUV and X-ray) is so efficient that it 
would destroy all PAHs (from regions in the disk where they could be excited) . As a possible path for 
PAH survival we suggest turbulent motions in the disk. They can replenish PAHs or remove them from 
the reach of hard photons. For standard disk models, where the surface density changes like and 
the mid plane temperature like r^^'^, the critical vertical velocity for PAH survival is proportional to 
r~^^* and equals ~5 m/s at 10 AU which is in the range of expected velocities in the surface layer. The 
uncertainty in the parameters is large enough to explain both detection and non-detection of PAHs. 
Our approximate treatment also takes into account the presence of gas which, at the top of the disk, 
is ionized and at lower levels neutral. 

Key words, dust, extinction - planetary systems: protoplanetary disks - infrared: stars - X-rays: stars 
- X-rays: ISM 



1. Introduction 

Infrarcxi emission bands of PAHs can be used as a 
probe of the UV environment. They are commonly 
seen in the ISM, but also in young stellar objects 
such as Herbig Ae/Be stars (Waelkens et al. 1996, 
Siebenmorgen et al. 2000, Meeus et al. 2001, Peelers 
et al. 2002, van Boekel et al. 2004). The observed 
emission can be explained in models of an irradi- 
ated disk (Habart et al. 2004, Visser et al., 2007, 
Dullemond et al. 2007a). 

ISO also looked at a few of the much fainter 
T Tauri stars but without a clear PAH detec- 
tion (Siebenmorgen et al. 2000). In the Evans et 
al. (2003) legacy program which employs the more 
sensitive Spitzer Space Telescope (SST), 3 out of 
38 T Tauri stars show PAH features (Geers et 
al. 2006). This corresponds to a detection rate of 
only 8% in contrast to almost 60% in Herbig Ac/Be 
stars (Acke & van den Ancker 2004). Similarly 
low rates for T Tauri stars are found by Furlan 
et al. (2006) who present 111 SST spectra in the 
Taurus-Auriga star forming region and speculate 
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that the absence of PAH resonances is due to the 
much weaker UV field compared to Herbig Ae/Be 
stars. Geers et al. (2009), on the other hand, argue 
that the PAHs arc simply under-abundant relative 
to the ISM. They also find that variations of the 
disk geometry, such as flaring or gaps, have only 
a small effect on the strength of the PAH bands. 
Clearing out gas and dust by planet formation in- 
side the disk could effectively remove PAHs. Indeed, 
inner gaps in disks are observed at radii between 
40-60 AU and at wavelengths between 20-1000/xm 
where the emission is dominated by large grains 
(Brown ct al. 2008, Geers et al. 2007b). However, 
in cases where PAH emission is resolved, it is ex- 
tended up to 15-60 AU, without sub-structure and 
inside the inner gap region (Geers et al. 2007a). The 
spatial extent of the PAH emission is also similar 
for T Tauri and Herbig Ae/Be stars. We therefore 
suggest that PAH removal by radiative destruction 
is dominant. 

Present radiative transfer models of the PAH 
emission from dusty disks consider only the stel- 
lar radiation field and no additional EUV or X-ray 
component (Habart et al. 2004, Geers et al. 2006, 
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Fig. 1. The spectral energy distribution of our T 
Tauri model star at 1 AU without foreground extinction 
(Eq. 1). The absolute luminosities of the components 
are given in Table 1. 



Visser et al. 2007, DuUemond et al. 2007). Their 
hard photons could, according to laboratory ex- 
periments (Ruhl et al. 1989, Leach et al. 1989a, b, 
Jochims et al. 1994) and theory (Omont 1986, 
Tielens 2005, RapacioU et al. 2006, Micelotta et al. 
2009), destroy PAHs. We discuss below their im- 
pact on the PAH abundance in the disks of T Tauri 
stars. 



2. Radiation components of T Tauri stars 

Our T Tauri model star has a total luminosity 
L^, ~ 2L0. Its radiation consists of a photospheric, 
a FUV, an EUV and an X-ray component. Their 
parameters are listed in Table 1 and are very simi- 
lar to those proposed by Gorti & HoUenbach (2008). 
The total spectrum is displayed in Fig.l. We point 
out that the FUV and EUV radiation are observa- 
tionally poorly constrained. 

The photosphere supplies most of the lumi- 
nosity whereas the FUV, EUV and X-ray radi- 
ation are much weaker and believed to originate 
from accretion onto the star and from chromo- 
spheric and coronal activity. The photosphere, the 
FUV and EUV component are approximated by 
blackbodies. We assume 4000 K for the photosphere 
and, following Stabler et al. (1980) and Calvet & 
GuUbring (1998), 15000 K for the FUV (pre-shock) 
and ^SxlO^K for the EUV emission (post-shock 
region) . 

The strength of the FUV and EUV radiation is 
determined by the accretion luminosity which we 
approximate by Lace = GM^M / R^. If i?, = 2R0 
and M* = 1 Mq are the radius and mass of the 
star, an accretion rate M = 10~^ yr~^ (Akeson 
et al. 2005) yields Lace ~ 0.01 L*. Higher values 



(Lacc/L* > 0.1), but with a large spread, are de- 
rived by MuzeroUe et al. (1998, 2003) from hydro- 
gen emission lines. However, as we show in sec- 
tion 5, such stronger fluxes have little influence on 
the stability analysis of PAHs. 

Preibisch et al. (2006) establish from Chandra 
observations (0.5-8 keV) a relation between the X- 
ray luminosity, L^, and the total luminosity con- 
firming the ROSAT results of Sterzik & Schmitt 
(1997). The ratio Lx/Li, is similar in rapidly ro- 
tating main-sequence stars and non-accreting T 
Tauri stars (^10^"^), but systematically lower by 
a factor in accreting T Tauri stars (Preibisch 
et al. 2006). Interestingly, in Herbig Ae/Be stars 
Lx/L^, is much smaller (^10~^ . . . 10~^, Stelzer et 
al. 2006) with values comparable to the Sun. The 
solar X-ray luminosity in the 0. 1-2.4 keV ROSAT 
passband lies during a solar cycle in the range 
10"^-^ < L^/L* < lO"^ '^ and is typical for G stars 
(Judge et al. 2003). 

X-ray fluxes are generally variable on timescales 
of hours to weeks and weaken during the evolu- 
tion of the T Tauri star. For example, half of the 
sources in the Taurus molecular cloud detected 
by XMM/Newton (0.3 - 7.8 keV) show variations, 
more at hard (> 0.5 keV) than at soft energies, and 
a quarter of them display flares (Stelzer et al. 2007), 
about once a week and lasting for a few hours. In a 
strong flare, more than 10'^^ erg are emitted and 
can reach 1% of the total luminosity. We assume 
up to 2 keV a power law spectrum oc (Giidel et 
al. 2007) and neglect harder radiation because the 
emission then steeply declines (oc v~^). 

Let Li be the frequency-integrated luminosity 
of the radiation component i (see Table 1) and Li^^ 
its spectral luminosity such that Li = J Li^^, dv. 
Dropping for convenience the index i, the flux (of 
component i ) at a distance r is 

where we included a screening factor e to ac- 
count for foreground absorption (by dust and gas). 
If K,u denotes the absorption cross section per car- 
bon atom, a PAH of Nc carbon atoms absorbs in 
one second (from component i) 
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Table 1. The four radiation components of our T Tauri model star. 
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(2) spectral shape (BB — blackbody), 

(3) mean energy of emitted photons, hV, 
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3. Cross sections 

As the light from the star enters the disk, it is atten- 
uated by gas and dust. The absorption cross section 
of gas depends on the ionization stage of the atoms 
which is determined by the balance between recom- 
bination and photo-ionization. By far the most im- 
portant atoms are, of course, hydrogen and helium 
with ionization potentials of 13.6 eV and 24.6 eV, 
respectively. Because the recombination rate is pro- 
portional to the square of the gas density which is 
high in the disk (section 5), the gas is ionized only 
in a thin surface layer {Ay < 0.001 mag, section 
2.6 of Gorti & Hollenbach 2008). We use atomic 
cross sections of Morrison & McCommon (1982) 
and Balucinska- Church & McCommon (1992) and 
solar element abundances. 

The dust cross sections are taken from the 
model of Kriigel ( 2006) which describes standard 
dust. For X-rays, the absorption efhciency calcu- 
lated from Mie theory must be corrected down- 
wards. Hard photons can eject electrons from the 
grain and as these carry away kinetic energy, only 
part of the photon energy is deposited in the dust 
particle. The threshold, Et, above which such a cor- 
rection is necessary depends on the grain size; de- 
tails are given in Dwek & Smith (1996). For a 10 A 
graphite particle, Et ~ 100 eV and the reduction 
factor is roughly proportional to 

The absorption coefficient of dust, K^^x, and of 
neutral gas plus dust, Kx = iiTgas.A + ^d,A) both 
per gram of disk material, are plotted in Fig. 2 for 
a dust-to-gas mass ratio of 1:130. Note that at the 
ionization threshold of hydrogen, /Tgas is almost 10^ 
times greater than K(\. 

With respect to the absorption cross section of 
PAHs, we assume = 7 x 10"^* cm^ per carbon 
atom when hi/ < 13.6 eV and scale at higher 
energies to follow the values of a graphite sphere of 
10 A radius (Dwek & Smith 1996). The maximum 
wavelength (in A) for PAH excitation is Amax,PAH = 
1630 + 370v^ (Schutte et al. 1993) resulting in a 
minimum photon energy of 2.3 eV for a PAH with 
A^c = 100 carbon atoms. 



3.1. PAH emission 

As the PAHs are transiently heated, their ex- 
citation is usually treated statistically. Following 
Guhathakurta & Draine (1989), let P{T) dT be the 
probability of finding in a large ensemble of PAHs in 
a steady state an arbitrary PAH in the temperature 
interval [T, T + dT] . The temperature distribution 
function P{T) is calculated in this method from 
a transition matrix {Afi). If denotes the PAH 
absorption cross section, the matrix element Afi re- 
ferring to dust heating from an initial enthalpy bin 
centered at Ui to a final one centered at Uj and of 
width MJf is, for a mono-chromatic flux, equal to 

t : else 

Examples of P{T) are displayed in Fig. 3 for 

mono-chromatic fluxes which cover almost the en- 
tire range encountered anywhere in the disk in 
terms of intensity, the flux ranges from F = 10 
to 10''ergs~^cm~^, as well as hardness, the pho- 
ton energy is between hu = 3.8 eV and IkeV. The 
approximate unattenuated fluxes of the four radia- 
tion components at a distance of 10 AU are listed 
in Table 1. Note that when F = n^hv is constant, 
the number of photons Uj decreases as the photon 
energy hi/ goes up. The power absorbed by one 
PAH is almost independent of the photon en- 
ergy as long as /ii/<100eV (see Fig. 4). By and 
large, when hv is fixed and F increases, the curves 
in Fig. 3 narrow and move to the right towards 
higher temperatures. When the radiation field is 
weak (F = 10 ergs^^cm"^) and the photons are 
soft {hv = 3.8 eV), the PAH absorbs about one pho- 
ton per day and there is plenty of time to cool down. 
In this case, the PAH virtually never exceeds the 
sublimation temperature Tg, (at which solid carbon 
gasify). For a wide pressure range (10^^ — lO^^dyn 
cm-2), Ts ~ 2000K for graphite (CRC Handbook 
of Chemistry & Physics 2005, Salpeter et al. 1977). 
In Fig. 3 we highlight the area where the tempera- 
ture is above Tg. For a PAH exposed to a strong 




Fig. 2. The mass extinction coefiicient per gram dislc material wlien tlie gas is neutral (Morrison & McCommon 
1983); the gas-to-dust mass ratio equals 130. 
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Fig. 4. The power W = Nc k, F absorbed by one PAH, 
with number of C atoms A^'c as indicated, in a mono- 
chromatic flux F = 10'' erg s~^ cm"'^ as a function of 
photon energy — Iw. 



radiation field {F = 10^ ergs~^cm~^), about a 
dozen of soft photons are absorbed within a coohng 
time and the temperature distribution function be- 
comes very narrow around 2000 K and sublimation 
is likely. For hard photons {hv > 50 eV) , the PAH 
undergoes, independent of the strength of the ra- 
diation field, extreme temperature excursions. This 
indicates that PAH become photo-unstable either 
by absorption of a single hard photon, or by soft 
photons if there are many of them. 



4. PAH destruction 

The abundance of PAHs is determined by the com- 
petition between formation and destruction pro- 
cesses under the specific environmental conditions. 
Underlying processes are discussed, for example, by 
Omont (1986), Voit (1992), or recently by Micelotta 
et al. (2009a). Here we only consider PAH de- 
struction by photons and generally assume that 
PAH formation is negligible. After photon absorp- 
tion, a highly vibrationally excited PAH may re- 
lax through emission of IR photons or, if suffi- 
ciently excited, lose atoms. The latter process is 
called unimolecular dissociation and is discussed 
for interstellar PAHs by AUamandola (1989), Leger 
et al. (1989), Le Page et al. (2003), Rapacioh et 
al. (2006), and Micelotta et al. (2009b). Laboratory 
studies of PAH dissociation which can be applied 
to astrophysical situations are rare (Jochims et al. 
1994). The photo-chemistry of PAHs is reviewed by 
Tielens (2005, 2008). 

4.1. Procedure 

In the disks of T Tau stars, the PAH abundance 
depends obviously on place and on time as the disk 
evolves. There is no general solution to the prob- 
lem and to extract numbers, we have to radically 
simplify it. We wish to find some estimate of the 
location where PAHs become stable against photo- 
destruction. To derive a procedure, we recall that 
although after absorption of an energetic photon its 
energy is immediately distributed over all available 
vibrational modes (AUamandola et al., 1989), the 
excitation of a particular atom fluctuates and occa- 
sionally it is pushed into the continuum and leaves 




temperature [K] 



Fig. 3. The temperature distribution P{T) of a PAH with 100 C atoms exposed to mono-chromatic radiation 
with hv — 3.8, 13.6, 50, 100, 300 eV and 1 keV. This set includes the mean photon energies of the four radiation 
components of the T Tauri star. The fluxes range from (top to bottom) _F = 10 to lO'^ erg s~^ cm~^. Shaded area 
marks temperatures above the sublimation temperature of graphite. 



the PAH. Quantitatively, the unimolecular dissocia- 
tion can be written in Arrhenius form. In a classical 
description, an atom of critical (Arrhenius) energy 
Eq detaches from a PAH of peak temperature Tp if 
the dissociation time 

idis - ^o"^e^"/^-^p (7) 

is shorter than the cooling time icooi- A charac- 
teristic value for the vibrational frequency is vq ^ 
10"'^'^ s~^. The "atom", which may besides H or C 



also be an atomic group like C2H2, needs the time 
^dis to overcome the critical internal barrier, Eq, 
which is similar but not identical to the chemical 
binding energy. Micelotta et al. (2009b) quote Eq 
of 3.2eV for H-loss, 4.2eV for C2H2, 7.5eV for pure 
C loss and 9.5eV for C2. For the ISM they find 
Eq = 4.6 eV and a somewhat larger value for PDR. 
The inverse of the dissociation time is the proba- 
bility that a certain atom leaves the PAH per unit 
time. 



The exponential term e^°/^'^<' in Eq.(7) in- 
creases very rapidly as T falls and meaningful val- 
ues (i.e. not too large ones) of idis are obtained 
only if Tp > 1500 K. Atoms will only detach when 
idis < icooi- A-s the cooling time at these tempera- 
tures is for astrophysical applications of order 1 s, 
independent of the PAH size, the dissociation cri- 
terion reads 



tdis ^ If 



(8) 



It leads to a minimum temperature for destruc- 



tion 



2dis = 



En 



fclnz^o 



(9) 



Assuming Eq ~ 5eV, one gets T^is — 2000 K. 
At this high temperature, the internal en- 
ergy of a PAH is reasonably well approxi- 
mated by 3NckTdis when also taking the pres- 
ence of H— modes into account. The minimum 
temperature T^is is related to a minimum energy in- 
put AE. The PAH is therefore unstable to photons 
with 

AE > hu^ = mckTii, = -^N^Eo ^ 0.1Arei;o(10) 

\nuo 

or when the number of carbon atoms 
2AE 



Nr. < 



[eV] 



(11) 



Micelotta et al. (2009b) find that a PAH with 
Nc = 50 requires an internal energy of about AE = 
24 eV to dissociate, which agrees with the above 
estimate (Eq. 11). 

The minimum energy input required for disso- 
ciation can either be delivered by absorption of i) 
many soft photons, with a total energy -Eabs ^ AE 
(Eq. 3), or ii) by a single hard photon, with energy 
hv > AE. If a photon heats the PAH to a peak 
temperature much above Tdis, more than one atom 
will detach. The first expulsion occurs momentarily 
(idis ^ Is). It consumes the energy Eq plus some 
kinetic energy i?kin for the liberated atom. The new 
PAH temperature follows from 

hi^- Eo- E^in = 3iNc-l)kT (12) 

This happens x times until T has dropped to Tdis, 

hi^ - x{Eo + Ei,in) = 3(iVe - x)kTai, (13) 

With Ekin ^ 0.5 eV, we estimate that the total 
number of freed atoms is 

^ Eo + i^kin - 3fcTdis ~ 5 [eV] 10 ^ ' 

For Nc = 100, an average EUV photon ejects nine 
atoms and an X-ray photon destroys the whole PAH 
(column (5) in Table 2). 



4.2. Disruption by Coulomb forces 

For completeness, we also mention disruption of 
PAHs by Coulomb forces. Double or multiple ion- 
ization of a PAH loosens the binding of the pe- 
ripheral H atoms as well as of the skeleton of car- 
bon atoms. The ejection of K-shell electrons by X- 
ray photons {hv > 284 eV) in combination with 
Auger electrons will amplify the process. Coulomb 
explosion is relevant mainly for small PAHs and ne- 
glected here. 



5. Conditions for PAH survival 

According to Eq. (10), PAHs are destroyed if the 
source emits photons of energy hv > O.INcEq, ir- 
respective of the distance to the star or its lumi- 
nosity. For Nc = 100, the critical photon energy is 
only 50 eV (Eq.lO). As T Tauri stars (or their jets) 
also radiate at X-rays and in the EUV, the surface 
of the disk should be devoid of PAHs unless a) the 
period over which hard photons are emitted is too 
short to destroy all PAHs; b) the PAHs are by verti- 
cal motions removed from the hard radiation before 
they are destroyed and there is an influx of PAHs 
from below; c) PAH destruction is compensated by 
PAH formation in the surface layer. The last ef- 
fect should, in a hard photon environment where 
PAHs and carbon atoms are ionized, be prohibited 
by Coulomb repulsion (Voit, 1992). 

5.1. Destruction time 

The above PAH survival condition under a) can 

easily be dismissed. To estimate the time for PAH 
removal, trem, by the radiation component i, we 
note that most of the radiation is absorbed on the 
disk surface in a sheet of vertical optical depth tj_ 
equal to the grazing angle ctg of the incident light. 
We call this sheet the extinction layer (of radiation 
component i) and denote its geometrical thickness 
ii (see Fig.5). To first order, the PAHs in the ex- 
tinction layer receive the stellar flux of Eq.(l) with 
~ 1. If the instability criterion of Eq.(lO) is ful- 
filled, trcm follows from xNjtrem = Nc, where X is 



(15) 



from Eq.(14), therefore 
Nc 

^rem — ^abs 
X 

With tabs from Table 2, one sees that even at 
100 AU, trcm is short compared to the duration of 
the T Tauri phase (-lO*^ yr, Bertout et al. 2007, 
Cieza et al. 2007). 

5.2. Exposure time and vertical mixing 

Next we consider the possibility that vertical mo- 
tions in the disk lead to a continuous exchange be- 
tween matter in the extinction layers, where almost 
all photons are absorbed and PAHs destroyed, and 




Fig. 5. Of each radiation component, ^^90% is absorbed in wliat we call the extinction layer. The optical depth 
from its bottom to the star is one and in vertical direction equal to the grazing angle Qg. The height of its lower 
boundary, zo, declines with radius, but its geometrical thickness is rather constant {£ ~ 0.5 H, see Fig.6). Vertical 
motions may replenish PAHs from below. 



the layers below where PAHs are shielded and dam- 
aged ones possibly rebuilt (Fig. 5). We assume that 
gas and dust are perfectly mixed in a mass ratio 
130:1. 

In a Keplerian disk that is isothermal in z- 
direction and in hydrostatic equilibrium, the gas 
density changes like 



p{z) = 



2 S 
H 



(16) 



Here S(r) is the surface density at radius r which 
is assumed to follow a power law, 



S(r) = 

and 
H{r) 



z) dz 



AU 



(17) 



m 



(18) 

is the scale height, ~ IMq the stellar mass 
and m the mass of a gas molecule. For the sur- 
face density, reasonable numbers are 7=1 and 
S = 200gcm~^ (Hartmann et al. 1998, Kitamura 
et al. 2002, Dullemond et al. 2002, Rafikov & Colle 
2006, Gorti & HoUcnbach 2008), although the var- 
ious estimates show considerable scatter. 

For the radial variation of the gas temperature 
in the opaque mid plane, T(r), we also adopt a 
power law, 

1-/3 



T{r) = To 



AU 



The mid plane is roughly isothermal in z be- 
cause the optical depth is high and the net flux 
zero. It is much colder than the extinction layers 
because it is not exposed to direct stellar heating. 
The radiative transfer in the disk, including the en- 
ergy equation, can be solved to any desired accu- 
racy even when the disk is very opaque (see section 
11.3.2 of Kriigel 2006). As long as the dust in the 
mid plane is optically thick to its own emission, 
the results for T(r) can be well approximated by 
putting in Eq.(19) = 0.5 and To 130 K (as also 
suggested by Dullemond et al. 2007b or Chiang & 
Goldreich 1997). 

Each extinction layer extends vertically from 
some value zq upwards to infinity (Fig. 5). We give 
it a finite thickness £ by demanding that, say, 90% 
of the photons are absorbed between zq and zq + i. 
If v± denotes the typical vertical velocity, for ex- 
ample, as a result of turbulence, PAHs are exposed 
to radiation for a time 



v± 



(20) 



This is also the mean residence time of a PAH 
in the extinction layer. For PAHs to survive, texp 
must be smaller than tiem Eq.(15). The height zq 
follows from 



(19) K / p{z) dz = Q!g 



(21) 




Fig. 6. The height, zq, of the bottom of the extinction layer and its thickness £ for the four radiation components 
(see Fig. 5 and Eq.(21), (22)) for a grazing angle ag = 3°. Due to the high gas densities, the EUV extinction layer 
is practically coincident with the ionization front and lies above the extinction layer of the other components. 
This implies that EUV radiation is absorbed first whereas the other components penetrate deeper. Neglecting 
ionization by the FUV component, the gas below the zo-\ine of the EUV radiation (top) is neutral. 



and i may be estimated from the condition that 
only 10% of the flux is absorbed above zq + i, 

/•oc 

K / p{z)dz = 0.1 ag (22) 

K is the mass absorption coefficient of gas and 
dust at the characteristic frequency of the particu- 
lar radiation component (see Table 2). Because p{z) 
changes rapidly, zq is rather insensitive both to ag 
as well as K. For Ug/K = 10~* . . . 10~^, one obtains 
Zq = 2.67?. . .5.7H. So in the V band, where ab- 
sorption is only by dust {K ~ 200 cm^ g~^) and for 
a grazing angle ag = 3°, one gets zq ~ {4. . .5) H 
and fi — t/H ~ 0.5. The height zq where an ex- 
tinction layer begins and its thickness are shown in 
Fig. 6; t is for all radiation components very similar 
{£, - H/2, i = 1,...,4). EUV photons are ab- 
sorbed highest up, their extinction layer lies about 
one scale height above the others. X-rays penetrate 
slightly deeper than photospheric or FUV photons. 

From Eq.(4), (15) and (20), one finds that ver- 
tical motions safeguard PAHs against destruction 
if 

OC 

with fc from Eq.( 10). When the scale height H 
is given by Eq.(18), Vcr oc r-^^+i'/^ = ^-3/4^ rp^^ 
values of t;cr, at a distance of 10 AU, are listed in 
Table 2 together with other quantities relevant to 
PAH survival. 



One expects the disk also to be turbulent. 
Turbulence may be driven by various processes such 
as shear flows in the disk (Lin & Bodenheimer 
1982), magneto-rotational instabilities (Balbus & 
Harley 1991) or velocity discontinuities at places 
where infalling matter (Cassen & Mossman 1981) 
or outflows (Elmegreen 1978) strike the disk sur- 
face. Until now it is not clear which type of tur- 
bulence dominates. We assume that for the size 
of the largest Eddies, ^cd, the average turbulent 
velocity, vt, grows linearly with the sound speed 
Cg. Various hydrodynamical 3-dimensional calcula- 
tions (Cabot 1996, Boss 2004, Johansen & Klahr 
2005, Fromang & Papaloizou 2006) support this 
view. The favoured parametrisation is vt = a'^Cs 
with q = 0.5. This choice has consequences on 
the Eddy scale £cd — a^~^H and the turn over 
time ted — ^cd/vt — a^~^''/r2K, with Keppler 
frequency fix = ^jGNh/R"^ (e.g. DuUemond & 
Domink 2004). Weidenschilling & Cuzzi (1993) use 
g = 1 so that the Eddy scale is about the pres- 
sure scale height, ^cd = H, and larger than the 
thickness of the extinction layer, fed > i ^ H/2. 
Estimated values for a are in the range from 0.0001 
up to 0.1 (DuUemond & Dominik 2004, Schrapler 
& Henning 2004, Youdin & Lithwick 2007). Taking 
a = 0.01 the Eddy scale is 10 times smaller for 
q = 0.5 than for g = 1, and in addition larger tur- 
bulent velocities are obtained with vt = \facs , sup- 
porting a faster transport of the PAH. Identifying 

in Eq.(23) with the turbulent velocity vt and as- 
suming a temperature dependence as in Eq.(19), we 
plot in Fig. 7 the vertical velocity as a function 
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Fig. 7. The critical vertical velocity for PAH survival Vc-c after Eq. (23) with respect to X-rays as a function 
of distance (full line). We identify the vertical velocity v±_ with the turbulent velocity vt. PAHs survive when 
v± = Vt > Vcr, else they are destroyed (shaded areas) by expulsion of atoms. The hatched strip between the 
dashed lines shows the range where nx is between 0.001 Cs and 0.01 Cs (where Cs is the sound velocity). 



Table 2. Quantities relevant to PAH survival. 





(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


component 


hV 








X 




l/H 


^abs 






eV 


lO-'^^A^ 














m/s 


photosphere 


2.7 


700 


< 1 


1.2 




4.4 


0.5 


5 s 




FUV 


4.2 


700 


< 1 


1.5 




4.5 


0.5 


120 s 




EUV 


97 


36 


~0.5 


110 


9 


5.3 


0.4 


9 days 


2900 


X-ray 


1100 


0.17 


< 1 


0.55 


100 


4.2 


0.5 


290 yr 


3.4 




(1) The mean energy of destructive photons, hV, from Eq.(5) at the bottom of the extinction layer (see Fig. 5); 
for photospheric and FUV photons we put Vc =Q and for EUV and X-rays we integrate for v > Vc (Eq. 10); 

(2) absorption cross section, k, per C atom at frequency V; 

(3) approximate mean degree of ionization of the gas in the extinction layer; 

(4) extinction cross section of gas and dust at frequency V normalized to Ky = 200 cm^ g~^; 

(5) number of expelled atoms, x, per absorption event from Eq.(14); 

(6) altitude of the bottom of the extinction layer in units of the scale height H; 

(7) thickness of extinction layer; 

(8) mean time tabs from Eq.(4) in which one photon is absorbed; 

(9) critical vertical velocity for PAH survival. 
Note: zo,^, iabs and Hcr refer to r = 10 AU. 



of radius. The figure also shows the critical velocity 
for PAH survival, v^, with respect to X-rays. Note 
that the critical velocity is for the X-ray radiation 
component insensitive to the particular choice of 
£^0 and (Eq. 10, 23). When v±_ = 0.01 Cs, PAHs 
can survive at distances r > 10 AU; when v±_ is 
considerably smaller than 0.01 Cg, they cannot. 

Critical velocities for PAH survival are much 
higher for EUV than for X-ray photons (Table 2). 
EUV radiation will therefore always destroy PAHs 
but, as depicted in Fig. 6, the EUV extinction layer 



is the topmost and below it, PAHs may survive and 
be excited. 

If PAHs are removed from the extinction layer 
before they are destroyed, they must, in order to be 
detected, at the same rate be injected from below. 
Therefore, the critical velocity can alternatively be 
expressed through 

/ p(z) dz ~ p(zo) Wcr (24) 

trom J zq 

which leads to similar values. We note that the mass 
reservoir below the extinction layer is sufhcient to 




Fig. 8. The energy Sabs (Eq.3), which is absorbed by a PAH of Nc — 100 carbon atoms, as a function of distance 
from the star. The PAH is exposed to the photosphcric and FUV radiation component described in Table 1. The 
dashed line refers to the top of the extinction layer (r — 0) and the full line to its bottom (t — 1). For Arrhcnius 
energy of Eq = 5, 7 and 10 eV, the minimum energy input AE (Eq.lO) for PAH dissociation is indicated by the 
dotted lines. 



sustain, over the lifetime of the disk tufe, the re- 
quired mass influx pVcr- 

5.3. PAH dissociation by soft versus hard photons 

The energy _Eabs (Eq.3), which is absorbed by a 
PAH of Nc = 100 carbon atoms, is shown in Fig. 8 
as a fimction of distance from the star. The PAH 
is exposed to the photospheric and FUV radiation 
component described in Table 1 and results are 
shown for the top (t = 0) and the bottom (t = 1) 
of the extinction layer. The minimum energy in- 
put AE (Eq.lO) for PAH dissociation depends on 
the choice of the Arrhenius energy Eq and is indi- 
cated for Eq = 5, 7 and 10 eV, respectively. In this 
picture, for Eq = 5cV and at the bottom of the 
extinction layer, PAHs are dissociated by soft pho- 
tons up to 1 AU. For X-rays, however, we find that 
PAH destruction occurs typically at distances up 
to ~ 10 AU or even larger (Fig. 7). Dissociation of 
PAH acts for soft (photospheric and FUV) photons 
on much shorter distances than for hard photons 
(X-ray component). 

6. Large grains 

Observations of T Tauri stars at millimeter wave- 
lengths (Testi et al. 2003, Lommen et al. 2007) 
and in the mid infrared (van Boekcl et al. 2003, 
Przygodda et al., 2003, Kessler-Silacci et al. 2006, 
Bouwman et al. 2008, Watson et al. 2009) suggest 
that grains in T Tauri disks are at least 10 times 
larger than those in the ISM. As such large grains 



may also be present in the top disk layer we esti- 
mate how this would affect the stability analysis of 
PAHs. We first note that in case of homogeneous 
mixing an increase in particle size would not alter 
the dust-to-gas mass ratio. 

Should the grains be much larger than the wave- 
length, the absorption coefficient per gram of dust, 
Kd.\, would decrease roughly like one over grain 
radius whereas the ratio K^^x/K^y would still 
roughly be given by the values in Table 2. For hard 
X-rays, on the other hand, K^^x is not sensitive to 
grain size. 

Therefore, if disk grains are on average ten times 
bigger and thus 10"^ more massive than interstel- 
lar ones, we expect that the height zq to which 
the stellar radiation components can penetrate (see 
Fig. 6) stays the same for X-rays, but also for EUV 
radiation because EUV absorption is due to gas, 
not dust. However, optical and FUV photons will 
reach farther down, about half a scale height, so 
that there may be a thin disk layer (~i?/4) where 
PAHs are shielded from X-rays and EUV photo- 
destruction and excited by optical or FUV radia- 
tion. 

7. Conclusion 

In the search for an explanation why most T Tauri 
stars do not exhibit PAH features and only a few do, 
we investigate which processes can remove PAHs 
from the surface layer of T Tauri disks and under 
which conditions they should be present. Clearing 
of PAH through interaction with planets seems not 



an efBcient process (Geers et al 2007b) and we show 

that PAH under abundance can be caused by ra- 
diative destruction. We use a fiducial model for the 
photon emission of the T Tauri star that includes, 
beside the photosphere, FUV and EUV radiation 
and an X-ray component. 

1. We introduce for each stellar radiation compo- 
nent the notion of extinction layer as the place 
where ^90% of the photons are absorbed. EUV 
photons are mainly absorbed by gas. X-rays by 
gas and dust alike and the photospheric and 
FUV component are only attenuated by dust. 
The extinction layer of all four components have 
a similar geometrical thickness and their bottom 
is at similar elevation zq, except for the EUV ex- 
tinction layer which lies higher up (Fig. 6). 

2. PAH may be radiatively destroyed, by unimolec- 
ular dissociation, where one or several atoms are 
expelled after photon absorption. 

3. Destruction by the photospheric and FUV radi- 
ation component (soft photons), increases with 
the strength of the radiation field and is very 
efficient below 1 AU. 

4. Hard photons can dissociate PAHs at all dis- 
tances and their efficiency grows with the hard- 
ness of the photons. Without some counter pro- 
cess, all PAHs (in layers where they can be ex- 
cited) would be destroyed within a time short 
compared to the lifetime of the disk. 

5. Although grains in the disk surface are presum- 
ably larger than interstellar ones, the stability 
analysis of PAHs would not change significantly. 

6. Therefore, in disks where PAHs are detected, 
there must be some survival channel. Because 
creation of PAHs in the extinction layer is too 
slow to compete with PAH destruction (Voit 
1992), we suggest vertical mixing as a result 
of turbulence. It can replenish PAHs or remove 
them from the reach of hard photons. 

7. For standard disk models, the minimum velocity 
for PAH survival is proportional to r^^^^ and 
equals '~5m/s at 10 AU. If turbulent velocities 
are proportional to the sound speed a velocity 
Vt > 5 m/s would imply vt/cs <; 0.01 as PAH 
survival condition. Theoretical predictions for 
this ratio have a large spread but in accordance 
with the observational fact that PAH features 
are usually absent it seems that generally the 
condition is not fulfilled. 

8. A higher PAH detection rate is found in Herbig 
Ae/Be stars. In our picture this is explained 
as their destructive hard radiation component 
is relatively weak {L^/L^ < 10"'', Preibisch et 
al. 2006) and because the intensity of the PAH 
emission from large distance from the star is 
larger given their higher optical luminosities. 
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